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Abstract 

We discuss the synchrotron radiation flux from the Galactic center in unsta- 
ble dark matter scenario. Motivated by the anomalous excess of the positron 
fraction recently reported by the PAMELA collaboration, we consider the case 
that the dark matter particle is unstable (and long-lived), and that energetic 
electron and positron are produced by the decay of dark matter. Then, the 
emitted becomes the source of the synchrotron radiation. We calculate 
the synchrotron radiation flux for models of decaying dark matter, which can 
explain the PAMELA positron excess. Taking the lifetime of the dark matter 
of 0(10^^ sec), which is the suggested value to explain the PAMELA anomaly, 
the synchrotron radiation flux is found to be 0(1 kJy/str) or smaller, depend- 
ing on the particle-physics and cosmological parameters. 
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1 Introduction 



Nowadays there is no doubt about the existence of dark matter in our universe. 
However, nobody knows what dark matter is. The answer for the question will 
give a great impact on various fields of physics including astrophysics, cosmology 
and particle physics. For astrophysics and cosmology, dark matter is the dominant 
source of gravitation, and it governs the structure formation of our universe. For 
particle physics, we cannot find the candidate for dark matter in the standard model, 
so that dark matter is now regarded as the first evidence of new particles beyond 
the standard model. 

Experimentally, many attempts have been made to understand the property of 
dark matter. Among them, the PAMELA collaboration has recently reported quite 
an interesting result; the anomalous excess of positrons with the energy of 20 — 100 
GeV in the cosmic ray [Ij. The excess is hardly understood by the conventional 
mechanism [21 [3] , namely the secondary production of positrons due to the collision 
between primary protons and interstellar medium in our Galaxy. 

Theoretically, a lot of scenarios for dark matter have been proposed to explain 
the anomalous excess of positrons. These scenarios are categorized into two types; 
positron productions from the annihilation of dark matter [H El [6l [71 El [3 |10l [11] 
and from the decay of dark matter [121 ilSl US HSl UHl [HI- In most of the scenarios 
in the former case, people have faced to the difficulty to explain the excess; a huge 
enhancement factor, called boost factor, to the production rate of positrons is needed, 
where the factor is from the inhomogeneity of the dark matter distribution in the 



vicinity of the solar system]^ On the contrary, in the latter case, the positron flux 
to explain the excess can be naturally obtained with the lifetime of dark matter 
much longer than the age of the universe. One might worry about the fine-tuning 
problem to have a long lifetime of dark matter, however, there are some scenarios to 
realize such a long lifetime, as we discuss in the following. 

Synchrotron radiation from the annihilation or decay of dark matter may give 
constrains directly to scenarios aiming to the explanation of the PAMELA positron 
excess. Since dark matter annihilates or decays into energetic positron (as well 
as electron) under the circumstance of magnetic fields in our galaxy, synchrotron 



*^There also be some exceptions, for instance, a huge annihilation cross section can be obtained 
if the Sommerfelt enhancement to the cross section is considered [HI [191 EQ] . A large boost factor 
is not mandatory in those scenarios 
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radiation is inevitably induced. Importantly, the WMAP collaboration has observed 
the radiation in the whole sky, so that the observation gives constraints on the 
scenarios of the electron and positron production due to the annihilation or decay 
of dark matter in the Galactic halo. The synchrotron radiation flux for the case of 
dark matter annihilation has been considered in past works [2T], • The purpose 
of this article is to study the synchrotron radiation from the decay of dark matter, 
and discuss implication to the scenarios of decaying dark matter which account for 
the PAMELA result. 

The organization of this paper is as follows. In the next section, we explain the 
cosmological scenario we consider. The formalism to calculate the synchrotron flux 
from the decay of dark matter is given in Section [31 Compared to the previous 
works, several improvements are realized. Then, we numerically calculate the flux 
for several dark matter models; the results are shown in Section m The summary of 
our study is given in Section O 

2 The Scenario 

In this section, we summarize the cosmological scenario we consider. As we men- 
tioned in the introduction, taking the anomalous excess of the positron fraction 
measured by the PAMELA [IJ seriously, we consider unstable dark matter scenario. 
Indeed, if the decay of dark matter produces energetic positron, it becomes the source 
of the cosmic ray positron, and may explain the excess of the positron fraction. 

In the following, we consider the case that the dark matter particle (which we 
denote X) is unstable, and that its decay produces energetic positron (as well as 
electron). We also assume that the energetic positron produced by the decay is the 
origin of the excess of the positron fraction observed by the PAMELA. So far, many 
possibilities of such an unstable dark matter has been discussed [T2l[T3l[TH [T5 |[T6l[T7] . 
The energy spectrum of the produced depends on the model; in some cases, only 
the monochromatic electron and positron are produced by the decay and, in other 
cases, the energy distribution of the is non-monochromatic. 

The positron flux from the dark matter decay is inversely proportional to the 
lifetime of dark matter tdm- The best-fit value of tdm to explain the PAMELA 
anomaly depends on the energy spectrum of the positron from the dark matter 
decay. For example, if the dark matter mostly decays as X — > e~^e~ , the PAMELA 
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anomaly is well explained when [13] 

r™^2.2xl0-secx(^^)-'. (2.1) 

For example, if the sneutrino lightest superparticle (LSP) in the supersymmetric 
model is dark matter and also if it decays via an i?-parity breaking superpotential 



interaction, this is the case 
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On the contrary, if dark matter is the unstable 
gravitino which dominantly decays into ifW^ final state (with if being charged 
lepton in i-th generation) [25l[26], is also produced by the decay of W^. In such 
a case, the best-fit lifetime becomes 

.™^8,5xl0-.ecx(^)-', ,2.2) 

In any case, in order to explain the PAMELA anomaly, the lifetime of O(10^^~^^ sec) 
is needed (when m-DM ~ 0(100 GeV)). In the following, we concentrate on the case 
that tdm is of this order. In addition, as we will discuss, taking these values of the 
lifetime and m-DM ~ 1 TeV, the (e~ + e+) flux becomes consistent with those reported 
by the ATIC [27] and the PPB-BETS |28] experiments, which have observed excess 
of the flux at the energy of a few hundred GeV. 

Once the electron and positron are emitted by the decay of dark matter in our 
Galaxy, they cause the synchrotron radiation and the emitted radiation may be 
observed. In the following, we calculate the flux of the synchrotron radiation from 
the central region of our Galaxy. 



3 Synchrotron Radiation Flux 

In this section, we formulate the calculation of the radio flux from synchrotron 
radiation generated by high energy electron (and positron). First, we discuss the 
propagation of the electron and positron produced by the decay in the Galaxy. 
Then, we summarize the formula for the calculation of the synchrotron radiation. 



3.1 Propagation model for the electron and positron 

The synchrotron radiation flux depends on the energy spectra of the electron and 
positron (at each point of the Galaxy) produced by the dark matter decay. We first 
discuss the propagation of the electron and positron in the Galaxy to understand 

*^Evcn if the sneutrino is dark matter, the direct-detection constraint can be avoided. See [231124] . 
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Model 6 Ko (kpcVMyr) L(kpc) i?(kpc) 



Ml 



0.46 



0.0765 



15 



20 



MED 



0.70 



0.0112 



4 



20 



M2 



0.55 



0.0060 



1 



20 



Table 1: Parameters for the propagation models of electron, which will be used in 
our numerical calculations. 

how the spectra are determined. For the electron and positron from the decay of 
dark matter, their spectra are expected to be the same, so we consider the electron 
spectrum. 

When the energetic electron is emitted, its trajectory is twisted in the so-called 
diffusion zone where the magnetic field is non- negligible. Because the magnetic field 
is entangled, the motion of the electron is expected to be well approximated by a 
random walk. Then, the propagation of the electron is described by the following 
differential equation: 



where fe is the electron spectrum (i.e., number density of the electron per unit 
energy), K{E) is the diffusion coefficient, b{E, x) is the energy loss rate, and Q{E, x) 
is the electron source term. 

Diffusion of injected electron is caused by the entangled magnetic field in the 
Galaxy. The function K{E) can be determined so that the cosmic-ray Boron to 
Carbon ratio and sub-Fe to Fe ratio are reproduced. In our analysis, we adopt the 
propagation model given in j29j, where the shape of the diffusion zone is approx- 
imated by a cylinder (with the radius R and the half-height L) and the function 
K{E) is parametrized as 



Parameters of three representative models, MED, Ml, and M2 models, are shown 
in Table [H The MED model is the best-fit to the Boron-to-Carbon ratio analysis, 
while the maximal and minimal positron fractions for E > 10 GeV are expected to 
be estimated with Ml and M2 models, respectively. 

The energy loss of electron is via synchrotron radiation under the magnetic field 
and inverse Compton scatterings with cosmic microwave background (CMB) and 



K{E)V^UE, x) + — [b{E, x)ME, x)] + Q{E, x) = 



(3.1) 
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infrared gamma ray from stars. Then, b{E, x) is given by the sum of the synchrotron 
energy loss rate Psynch and the inverse Compton energy loss rate Pic- 

b{E, X) = Psynch + Pic = Psynch [l + '^IC/synch (^)] • (3.2) 

For simplicity, we neglect the position dependence of the magnetic flux density B 
in the diffusion zone; then, Psynch becomes uniform in the diffusion zone. Since the 
magnetic flux in the Galaxy is not well understood, we take several values of B to 
see how the synchrotron radiation flux depends on the magnetic field in the following 
numerical analysis. The synchrotron energy loss rate Psynch is given by [30l 



synch 



^-^-^ ~ 3.4 X 10-^^ GeVsec-i x ( ' ( ' , (3.3) 



where e is the electron electric charge, eo is the permittivity of free space, rrie is the 
electron mass, and c is the speed of light. In addition, rjc/synch is expressed as 

2 f/rad(^) „ „ / f/rad \ f P ^ 



no,.,„o„(J) = 3(^57^-2.7x^j^^-— J , (3.4) 

where t/rad is the radiation density, and /io is the permeability of free space. In order 
to take into account the position dependence of U^ad, we use the results given in [32] 
where the radiation density on the Galactic disc (i.e., ?7rad('"|j, = 0)) and that on 
the surface of ry = 4 kpc (i.e., Uj-ad{f\\ = 4 kpc, z)) are given, where z and are the 
height and the radial distance in the cylindrical polar coordinate. We assume the 
following formula for the radiation density (unless otherwise mentioned): 

= 777TT^f^'^(^ll'0) + f^CMB, (3.5) 
f4(4kpc,0) 

with UcMB — 0.3 eVcm~^ being the CMB contribution to the radiation density. In 
addition, f4 = ^rad — f^cMB is the radiation density from stars, which is determined 
from the results given in |32]. Then, one can see that f/rad varies as 0.3—10 eVcm~^ 
from the rim to the center of the Galaxy. Notice that, in many previous works, t/rad 
is approximated to be a constant. However, the position dependence of Uj-ad gives 
a significant effect on the resultant synchrotron radiation flux because the inverse 
Compton scattering may become the dominant energy-loss process. 

The electron source term is given by electron injection rate and dark matter 
distribution in the Milky Way Galaxy halo as 

Q(s,f) = ^£2MW:§, (3.6) 

tdm mnu dE 
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where pdm is energy density of dark matter and dNe/dE is energy distribution of 
electron from the decay of single X. For the dark matter distribution Pbm, we use 
the isothermal halo density profile: 



^2 I ,^2 

Phaio(r-) = P0 7', J ' (3-7) 

core ' 

where ~ 0.43 GeV/cm^ is the local halo density, rcore — 2.8 kpc is the core radius, 
Tq ~ 8.5 kpc is the distance between the Galactic center and the solar system, and 
r is the distance from the Galactic center. 

Before solving the diffusion equation (13. ip . it is instructive to estimate the typical 
propagation length of the electron per time scale of the energy loss, which we denote 
{1{E)). For the electron with energy E, the time scale of the energy loss is ~ E/b{E), 
and hence we obtain {1(E)) ~ a/ KE/b. For the propagation model given in Table 
[1], {1(E)) becomes a few kpc for ~ 1 GeV (where we have used the energy loss 
rate given in (13.21) ). and it becomes shorter as the energy increases. For such a scale, 
the change of the dark matter density is not significant (except for the central region 
of the Galaxy). In particular, if the effect of the propagation is negligible, we can 
omit the term proportional to K(E) in the diffusion equation and, in such a case, 
we obtain the following approximated formula for the electron spectrum: 

(■(local) / Z7 _ ^ Pdm(x) Ye(> E) 



fr^^^(E, X) = (3. 
^DM "^DM 0(E,X) 



where 



oo 



Thus, in such a case, the energy spectrum at the point x is almost determined by 
the local dark matter density at the same point. 

In Fig. [H we plot the electron flux $e = at the solar system (i.e., r = Tq and 
z = 0). Here, we consider the case that Br(X —>■ e^e~) = 1, taking moM = 200 GeV 
and 600 GeV, tdm = 5 x 10^^ sec, and B = 3 /iG; then, the energy distribution in 
Eq. dM]) is 

dN 

^ = 5(E-m^M/2). (3.10) 
dE 

In addition, in depicting Fig. [T], we use f/rad — 0.9 eV cm~^ for simplicity. One can 
see that the flux fj^^"'^^^^ -well agrees with the predictions of the MED and Ml models 
for E > 10 GeV, which confirms our earlier discussion. On the other hand, one 
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Figure 1: Electron fluxes in MED (dashed), Ml (dot-dashed), and M2 (dotted) 

models, and that from Jp"^*^') (solid), for the case that Br{X — e"^e"") = 1. Here, we 
take moM = 200 GeV (top) and 600 GeV (bottom), tdm = 5 x 10^6 sec, B = 3 /xG, 
and ?7rad = 0.9 eV cm"^. 



can observe a significant difference between /i ' and the fiux in the M2 model. 
This is due to the smallness of half-height of diffusion zone in the M2 model, which 
is 1 kpc; with such a small half-height, the electron can escape from the diffusion 
zone. Given the fact that the MED model is the best fit to the cosmic-ray Boron 
to Carbon ratio analysis, we adopt the MED model as our benchmark model of the 
electron propagation. 
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3.2 Synchrotron radiation: formalism 

Synchrotron radiation energy density per unit time and unit frequency is expressed 
as 



L^{x) = j dE V{u,E)fe{E,x). (3.11) 



Here, Vlu, E) is synchrotron radiation energy per unit time and unit frequency from 
single electron with energy E, which is given by 



1 \/^p^ B 

V{u,E) = -—^ F{u/u,), (3.12) 

Aireo nieC 

where v^. is critical frequency defined as 

3eE2S 

""^^^^^^ ^^-'^^ 

and 

POO 

F{x)=x dyK,/,iy), (3.14) 



with Kz being the modified Bessel function of 2;-th orderly The function F{v/v^ has 
a peak at z/ ~ 0.29z/c- Adopting the Galactic magnetic flux density of i? ~ 3 ;uG, we 
can see that the synchrotron radiation in the observed frequency band of the WMAP 



(i.e, 22 — 93 GHz) is from the electron with the energy of -E ~ 10 — 100 GeV 
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For 



the electron in such an energy range, as we have seen, the electron spectrum with the 
MED model, which is our benchmark model, is well approximated by f^^""^^^^ given 
in Eq. fl3.8p . Thus, we use in our numerical calculation. 

In order to calculate the observed radiation energy flux, we integrate L^{x) along 
the line of sight (l.o.s.), whose direction is parametrized by the parameters 9 and 



^•^This formula, as well as energy loss rate given in Eq. (|3.3p . is for the case that the velocity of 
electron is perpendicular to the magnetic field. Of course, such situation is not always satisfied in 
the Galaxy. Therefore, considering the fact that the injected electron is isotropic and the magnetic 
field is entangled, we take the B in the formula as mean magnetic flux density which the electron 
feels effectively. 

'^"'in the previous works [331 [Ml [351 [351 ISZ] , P{^) is approximated as 

F(x) ~ ^5{x-Q.2Q). 
9%/3 

However, we found that such an approximation is not good in particular when Ye(> E) is suppressed 
for E > y/ATTmlc^i'/3eB. 
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0, where 9 is the angle between the direction to the Galactic center and that of the 
line of sight, and (j) is the rotating angle around the direction to the Galactic center. 
(The Galactic plane corresponds to = and vr.) Then, the synchrotron radiation 
flux is given by 



J,(0,0) = ^ / dlL,{l). (3.15) 

'l.o.s. 
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We note here that we do not consider the synchrotron radiation from the very central 
region of the Galaxy and that the effect of the self absorption is unimportant |36] . 
Using Eqs. (13.81) and (I3.12p . we obtain the energy flux from synchrotron radiation as 

, r (3.16) 

327r^ moMTDM Jl.o.s. 1 + nc/synch(x) J Vc 

Notice that, adopting the approximation of the constant magnetic flux in the Galaxy, 
the line of sight and energy integrals f actor ize. 



4 Numerical Results 

Now, we show the synchrotron radiation flux. We numerically evaluate the flux given 
in Eq. (13.161) for several dark matter models. In the following, we show the results 
for = |, i.e., JuiP, 7r/2). The frequency of the radiation is taken to be some of the 
WMAP frequency bands: u = 22, 33, and 61 GHz. 



4.1 Leptonically decaying dark matter 

Let us start with the simplest case where the dark matter particle dominantly decays 
as X — > e^e~. This is the case if, for example, the sneutrino field in supersymmetric 
model is the LSP, and also if the LjLi^i-type i?-parity violating superpotential exists 
(where Lj and Ei are superfields for SU{2)l doublet- and singlet-lepton superfields, 
respectively) [13j. In this case, the energy distribution in Eq. (13. 6p is given by Eq. 
(13.101) . As we have mentioned, with such a flux, the PAMELA anomaly can be 
well explained if the lifetime of dark matter is properly chosen. (See Eq. (12. ip .) In 
addition, when is larger than ~ 100 GeV, the dark matter contribution to the 
electron and positron fluxes for the energy range of 100 GeV < -E < |mDM are 
approximately given by 

9 




Figure 2: Synchrotron radiation fluxes at z/ = 22 GHz as functions of dark matter 
mass for angle 9 = 5°, 10°, 15°, and 20°. Here, the dark matter is assumed to decay 
only into e^e~ pair. Here, we take = 5 x 10^^ sec, and show the cases of B =1, 
3, 10 fiG (from the bottom to the top) for each figure. 



where the MED model of the propagation is used here. (The flux for E > ^moM 
vanishes.) Then, if mDM ~ 1 TeV, a bump at -E ~ 500 GeV shows up in the (e' + e"*") 
flux, taking the best-fit lifetime to explain the PAMELA anomaly. The behavior is 
consistent with the excesses of the (e~ + e"*") flux observed by the ATIC [27] and the 
PPB-BETS [28] experiments. 

In Fig. [2], we plot the fluxes for u = 22 GHz as functions of mDM, for tdm = 
5 X 10^^ sec, which is of the same order of the lifetime suggested from the PAMELA 



anomaly!^ The angle 6 is taken to be 5°, 10°, 15°, and 20° (which give the minimum 
distance from the Galactic center to the line of sight of 0.7, 1.5, 2.3, and 3.1 kpc, 
respectively), and B = 1, 3, and 10 fiG. We also show the fluxes for u = 33 and 
61 GHz in Figs. [3] and HI respectively. 

As one can see, the synchrotron radiation flux from the decay is 0(1 kJy/str) or 
smaller in wide parameter region mDM ~ 100 GeV — 1 TeV. The flux J^, decreases 



*^The synchrotron radiation flux is inversely proportional to tdm (see Eq. (|3.16[) ). In order to 
obtain the flux for different value of tdm, one can simply rescale the flux given in the figure. 
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Figure 3: Same as Fig. [21 except for u = 33 GHz. 



as the frequency increases. In addition, J^, decreases as B decreases; this is because, 
as the magnetic field becomes weaker, the e^s lose more energy by inverse Compton 
process, resulting in a suppression of the synchrotron radiation flux. 

Radiation flux from Galactic center region has been observed by the WMAP 
for frequency bands of 22, 33, 41, 61, and 93 GHz [3HI ES]- Since then, intensive 
analysis has been performed to understand the origins of the radiation flux. (For 
recent studies, see [381 1211 [221 [39].) Most of the radiation flux is expected to be 
from astrophysical origins, such as thermal dust, spinning dust, ionized gas, and 
synchrotron radiation, which have been studied by the use of other survey data 
[ini [m [12]. With the three-year data, the WMAP collaboration claimed that the 
flux intensity can be explained by the known astrophysical origins |43j. On the 
contrary, Refs. [211 [221 [SB] also studied the WMAP three- year data, and claimed that 
there exists a remnant flux from unknown origin which might be non-astrophysical; 



the remnant flux is called the 'WMAP Haze'' 
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However, no clear indication of the 



[211 [22], as a possible explanation for the WMAP Haze, annihilation of dark matter is 
proposed. If energetic charged particles are emitted via the dark matter annihilation, they become 
another source of the synchrotron radiation. According to the study, the synchrotron radiation 
flux turns out to be comparable to the WMAP Haze with the annihilation cross section {av) ^ 
O(10~^^ cm^sec"^), which is consistent magnitude to explain dark matter abundance. 
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Figure 4: Same as Fig. [21 except for u = 61 GHz. 



WMAP Haze from unknown source was reported by the WMAP collaboration after 
five-year data [39] . 

The existence of the WMAP Haze seems still controversial, and the detailed 
studies of the WMAP Haze using the data is beyond the scope of our study. Here, 
we adopt the flux of the WMAP Haze suggested in f2T[ [22] as a reference value. 
Our numerical calculation shows that the predicted flux, which is 0(1 kJy/str), is 
comparable to the flux of the WMAP Haze given in [2n [22] . As we mentioned, 
since the the existence of the exotic radiation flux of this size is controversial, it is 
difficult to confirm or exclude the present scenario without better understandings of 
the sources of Galactic foreground emission. 

4.2 Gravitino dark matter 

Next, let us consider the case that the gravitino (which we denote ipfj^) is the LSP and 
is dark matter. Even if the gravitino is the LSP, it becomes unstable if the i?-parity 
is violated. In particular, if the effect of the i?-parity violation is small enough, the 
lifetime of the gravitino becomes longer than the present age of the universe and it 
is a viable candidate for dark matter [44j. Production of the gravitino may be due 
to the thermal scattering |45] , the decay of the lightest superparticle in the minimal 
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supersymmetric standard model sector |l6l STJ HHl HHl En], and the decay of the 
inflaton field [SI]. We do not specify any particular scenario, but just assume that 
the density parameter of the gravitino has somehow become consistent with the dark 
matter density. 

In this study, we consider the following i?-parity violating operators: 

£rpv = BiLiH^ + ml^^LiHl + h.c, (4.2) 

where Lj is left-handed slepton doublet in i-th generation, while and are up- 
and down-type Higgs boson doublets, respectively. (We work in the basis where the 
bi-linear i?-parity violating superpotential vanishes, which is realized by an appropri- 
ate redefinition of the Higgs and lepton-doublet chiral superfields.) With the above- 
mentioned i?-parity violating operators, the anomalous excesses of the positron [521 [T] 
and 7-ray [Mj in the cosmic ray can be simultaneously explained when the parame- 
ters Bi and m| are properly chosen]^ 

The fiavor of the final-state lepton depends on the sizes of the i?-parity violating 
coupling constants. In this analysis, for simplicity, we assume that the gravitino 
decays only into leptons in one of the three generations. Then, the gravitino decays 
as ipti ~^ ^/ ^f^^ ^/^' ^"^^ (where / = 1, 2, or 3 specifies the generation 
of the final-state lepton). In particular, the process — > I'fW^ has the largest 
branching ratio (when the process is kinematically allowed). Notably, from this 
process, the charged lepton is directly produced, which results in an energetic 
electron or positron. If = e^, such an has energy of ("^3/2 ~'"^h/)/2"^3/2, where 
m3/2 is the gravitino mass. In addition, the decays of the weak- and Higgs-bosons 
also produce e^. As we mentioned in Section [21 the predicted positron fraction well 
agrees with the PAMELA result if the lifetime of the gravitino is given by Eq. (12.21) . 
In addition, in such a case, when 7/13/2 is large enough, the electron and positron 

"^^It is often the case that, when the decay of dark matter produces hadronic objects, hke qq pair, 
the cosmic-ray anti-proton flux becomes too large to be consistent with observations. However, the 
anti-proton flux is sensitive to the model of cosmic-ray propagation 25]. In addition, the ratio of 
the number of the qq pair to that of the positron is ~ 6 in the present scenario, which is about 4 
times larger than the ratio for the scenario proposed in jl4j where the heavy gauge boson for an 
exotic U{1) gauge symmetry, which has a small kinetic mixing with the U{1)b-l gauge boson, is 
dark matter. The cosmic-ray anti-proton flux in such a scenario is order of magnitude smaller than 
the observed anti-proton flux. Thus, we expect that the constraint from the anti-proton flux does 
not exclude the present scenario. More detail of this subject will be discussed elsewhere [51] . 
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fluxes for the energy range of 100 GeV ^ E < ^m^M are approximately given by 

E^<l>,± ~ 3.0 X 10"^ GeV/cm^ sec str x ( ^^A- ^ f^^^) ^ , (4.3) 

' V8-5 X 1025 secy VlTeV/ ^ ^ 

where is the hfetime of the gravitino. Then, as in the previous case, there exists 
a bump in the (e~ + e"*") flux at ii^ ~ 500 GeV if we take ~ 1 TeV (and 

the PAMELA best-flt value of the lifetime). This may be the origin of the excess 
of the (e- + e+) flux recently observed by the ATIC [27] and the PPB-BETS [28] 
experiments. (For detail, see [S]-) On the contrary, if = /i^ and r^, from the 
decay of becomes non-monochromatic. 

The electron and positron produced by the gravitino decay become the sources of 
the synchrotron radiation. In our numerical calculation, we include all the relevant 



#8 



Then, 



interaction terms to calculate the partial decay rates; for detail, see p6] 
for a precise calculation of the energy spectrum, we use the PYTHIA package [55] 
to calculate dNf./dE. 

First, we show the results for the case where / = 1 (i.e., = e^), for which J^, 
is maximized for the gravitino dark matter case. In Fig. [5l we plot the fluxes for 
V = 22 GHz as functions of the gravitino mass, taking r3/2 = 5 x 10^^ sec. The angle 
9 is taken to be 5°, 10°, 15°, and 20°, and -B = 1, 3, and 10 /iG. One can see that 
the flux is maximized when the dark matter mass is slightly larger than 100 GeV. 
This is due to the fact that, when ttidm ~ 100 GeV, Vc for from the decay of X 
deflned in Eq. (13.131) becomes comparable to v ~ 0(10 GHz). 

The second case is that the primary lepton produced by the gravitino decay is in 
the second generation. The results are shown in Fig. O The synchrotron radiation 
flux becomes smaller in this case, compared to the case of / = 1. This is because 
the produced by the decay of //^ is less energetic than that directly produced by 
the process ip^ e^W^. We have also studied the case that = t^. As expected, 
the radiation flux becomes more suppressed in such a case. 

For the gravitino dark matter case, the synchrotron radiation flux is again of 
the order of ~ 1 kJy/str or smaller. As we have mentioned, we expect that the 
synchrotron radiation flux of this size is not excluded by the results of the presently 
available observations. The flux in the gravitino dark matter case is smaller than 
that for the dark matter which decays into e~^e~ pair (if the lifetime is flxed). This 



#^There are typos in Eqs. (3.18) and (3.19) of [26 . The factors in front of the terms proportional 

- not ^ 
3 ' ^^^^ 2 ■ 



to the function G should be |, not |. We thank L. Covi for pointing out these typos. 



14 



X 
3 



w 2 



X 
3 








100 100 1000 

Gravitino mass (GeV) Gravitino mass (GeV) 



Figure 5: Synchrotron radiation fluxes at z/ = 22 GHz as functions of gravitino mass 
for angle 6 = 5°, 10°, 15°, and 20°. The final-state lepton in the gravitino decay is in 
the first generation. Here, we take ts/2 = 5 x 10^^ sec, and show the cases of B =1, 
3, 10 fiG (from the bottom to the top) for each figure. 



is due to the fact that the total amount of energy carried away by is smaller in 
the case of the gravitino dark matter. 

5 Summary 

In this paper, we have discussed the synchrotron radiation flux from the Galac- 
tic center in unstable dark matter scenario. Motivated by the recently reported 
PAMELA positron excess, we consider unstable dark matter whose decay produces 
energetic electron and positron; if the lifetime of dark matter is O(10^^~^^ sec), the 
observed positron flux can be well explained in such a scenario. Once the energetic 

is produced in the Galaxy, it becomes the source of the synchrotron radiation. 
As discussed in Section [3l in our study, we have used a formalism with several im- 
provements to calculate the synchrotron flux. Then, we have numerically calculated 
the flux for models which can well explain the PAMELA anomaly. 

Assuming tdm ~ 0(10^^ sec) to explain the PAMELA anomaly, we found that 
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Figure 6: Same as Fig. [5l except for the case that the final-state lepton in the 
gravitino decay is in the second generation. 

the synchrotron radiation flux from the dark matter decay is expected to be Ji, ~ 
0(1 kJy/str) (or smaller). As we have mentioned, the existence of the exotic ra- 
diation flux of this size is controversial. However, it should be also noted that the 
synchrotron radiation flux is inversely proportional to the lifetime of dark matter. In 
order to explain the PAMELA anomaly, the smaller value of tdm is preferred as ttidm 
becomes lager. (See Eqs. (12. ip and (12. 2p .) Thus, if we consider too large mDM, it 
may become difficult to explain the PAMELA anomaly in the decaying dark matter 
scenario once we take into account the WMAP observation of the radiation from the 
central region of our Galaxy. 

In any case, without better understandings of the sources of Galactic foreground 
emission, it is difficult to confirm or exclude the scenario of explaining the PAMELA 
anomaly in the decaying dark matter scenario (with ~ 0(100 GeV)). A more 
detailed study of the Galactic emission may provide a significant test of the unstable 
dark matter scenario. Once the existence of the WMAP Haze will be somehow 
confirmed, decaying dark matter may be a viable candidate of its origin. 
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